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Abstract 

o 

^SJ ' The goal of this paper is to understand from a mathematical point of view the 

magnetic confinement of plasmas for fusion. Following Frenod and Sonnendriicker [7] , 
^ ■ we first use two-scale convergence tools to derive a gyrokinetic system for a plasma 

submitted to a large magnetic field with a slowly spatially varying intensity. We 
formally derive from this system a simplified bi-temperature fluid system. We then 
investigate the behaviour of the plasma in such a regime and we prove nonlinear 
stability or instability depending on which side of the tokamak we are looking at. In 
Oh I our analysis, we will also point out that there exists a temperature gradient threshold 

i beyond which one can expect stability, even in the "bad" side : this corresponds to the 

rH ' so-called H-mode. 

1 Introduction 

psj , 1.1 Magnetic confinement for plasmas 
> 

' Fusion is undoubtly one of the most promising research fields in order to find new sources 

<N . of energy. For the time being, magnetic confinement fusion represents one of the two main 

, approaches (the other one being inertial confinement fusion). The principle consists ba- 

I sically in using a magnetic field in order to confine the very high temperature plasma. 

Good confinement is absolutely compulsory since the plasma could otherwise damage the 

0^ . surrounding materials. 

o 

> 

•i-H 

X 




A first step towards confinement is to use a tokamals^, i.e. a torus-shaped box and 
consider a large purely toroidal magnetic field B, in other words 3 = ^6^^ with e > 
small. One can formally show that at leading order in e, particles oscillate around the 
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magnetic field lines. The drawback of this technique is that there are in fact many drifts 
appearing at higher order, some due to the geometry of B and one we are specifically 
concerned with, which is called the electric drift or E x B drift: 



VE 



EAB 



where E denotes the electric field. 

Since the electric field is induced by the plasma itself, one can not precisely predict its 
qualitative behaviour and thus this drift may prevent us from getting a good confinement 
property : if we wait long enough, particles may stop to perfectly turn around the torus 
and start drifting toward the edge of the tokamak. In order to overcome the effects of the 
electric drift, the idea is basically to take advantage of the other drifts due to the geometry 
of B. 

In the present paper, we make the assumption that the ions of the plasma are at ther- 
modynamic equilibrium and we describe the distribution of electrons by a kinetic equation. 
For the sake of simplicity, we restrict to the 2D problem in the plane orthogonal to B, in 
order to understand the behaviour of the particles in the slice. We take a magnetic field 
given by 

B 



B 



with e > a small parameter and B to be fixed later. We consider the Finite Larmor 
Radius scaling (see |7j for a reference in the mathematical literature) which consists in 
considering a characteristic spatial length with the same order as the Larmor radius (which 
is of order e). This scaling allows for a better description of the orthogonal motion and is 
expected to make the electric drift appear in the limit e — )• 0. The density /^(t, x,f) (with 



t > 0, X E T^, u G R^) of the electrons is then given by the following dimensionless Vlasov 



Poisson system : 



fe,\t=Q = /O 







(1.1) 



fedv - 1 



We denote x 
vector A - 



Xi 



Vl 



in the local orthogonal basis (see figure [T]). For any 



X2j \V2, 

we denote A-^ 

Following Grandgirard et al. ([9]), we consider the explicit formula for B : 

Rq Rq 



A2 



A2 

-AiJ- 



B 



Rq + er cos 9 Rq + exi 



(1.2) 



denoting by Rq the small radius of the torus. (We recall that the characteristic spatial 
length is of order e) 

We consider that i?o ~ 1 j consequently at first order in e we get: 



B = l-exi 



(1.3) 



2 




Figure 1 : A slice of tokamak 



leading to the following system: 



dtfe + 7.V,./e + (E, + 
fe,\t=0 = fo 



xiv^).V.,f, = 



(1.4) 



fdv - 1 



We will see that taking an inhomogeneous intensity for the magnetic field, even at order 1 
in e, leads to a quite different behaviour for the plasma. 

Indeed, in the limit e — )• 0, we can derive rigorously another kinetic system which is 
qualitatively close to the following one (see sections [2] and [3|) : 



f dtf -\\v\''d^,j + E^.v,f = Q 

f\t=o = fo 



E 



-v.v 



(1.5) 



fdv-1 



Observe here that E-^ corresponds to the electric drift Ex B that we mentioned earlier; the 
additional drift Vd = —\v'^e2 is due to the inhomogeneity of the magnetic field intensity. 
The remarkable point is that this drift has a fixed direction; it makes the particles "fall" 
toward the "bottom" of the slice. At this point of the modeling, we now have to distinguish 
between the plasma-core and the plasma edge (see figure 1), the only difference between 
the two we are concerned with, being that the core is much hotter than the edge. This 
means from a kinetic point of view that the velocities are much smaller in the edge. 

We now divide the slice into two areas: we denote the part j;i > the "bad curvature" 
side and the part xi < the "good curvature" side: indeed, we expect the plasma in the 
"good curvature" side to be well confined, while the plasma in the "bad curvature" region is 
badly confined. This behaviour can be easily predicted with the following heuristic study 
in the "bad curvature" side: 
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Particles in the hot plasma drift faster (left figure), so if there is any perturbation (right 
figure), there appears a separation of charge creating an electric field E, which entails a 
drift E-^ that accentuates the perturbation: in other words, the equilibrium is unstable. 
This discussion is part of the folklore in plasma physics for tokamaks and this instability 
is recognized to be one of the main sources of disruption for the plasma. 

In the other hand one can lead the same qualitative analysis in the "good curvature" 
side and show in this case stability. 



1.2 Objectives and results of this paper 



In this paper, following the previous heuristic argument, we will particularly focus on 
system (II. 6p . which is a kind of simplified "bi-fiuid" version of (II. 5p . We consider that the 
plasma is made of two mixable phases, one being the hot plasma (with constant temperature 
T"*" and density p'^{t,x)) and the other the cold plasma (with constant temperature T~ 
and density p~{t,x)) . Of course, hot and cold means that > T~ . 



dtp+ - T+d,,p+ + E^.V,p+ = 
dtp- - T'd^,p+ + E^.V^p- = 
E = -V^V 
-A^V = p+ + p- - I 
y = on xi = 0, L 

(p+,p")|t=o = (Pc[,Po) with J Po +Po 

for t > 0, X £ [0, L] X M/LZ and L is the size of the box. 
The temperature T(i, x) of the plasma is given by: 



(1.6) 



Tit,x) 



p+{t,x)T+ + p~{t,x)T- 
p+{t,x) + p-{t,x) 



(1.7) 



(We refer to Section [3] for more details.) 

Unfortunately, we were able to derive this system only formally from system (jl.4p 
and had to make some physical and mathematical approximations. These are precisely 
explained in Section [3l 

We observe that this system shares structural similarities with 2D Euler equations in 
vorticity form, which describe an incompressible inviscid fiuid. 

Plasma -f-)- Fluid 

density /? — 1 ^ vorticity oj 
(rotated) electric field E-^ = V^A"^(/9 - 1) o velocity u = V^A^^tJ 
Hence, our system can be seen somehow as a Euler system with two kinds of vorticities. 
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Such an analogy between strongly magnetized plasmas and bi-dimensional ideal fluids 
has been observed for a long time by physicists (for instance, see |14)). We mention that 
the convergence towards 2D Euler in strong magnetic fields regimes (but different from the 
one studied here) was rigorously established by Golse and Saint-Raymond [H] and Brenier 

m- 

Let us now define precisely the stability and instabihty notions that we will work on until 
the end of the paper. One should be aware that for such infinite-dimensional dynamical 
systems, the choice of the norm is particularly important. 

Definition. Let ^ he a solution to M.b]) . This solution is said to he stahle with respect to 
the X norm if for any r] > 0, there exists 6 > such that: for any solution p to U-6\) . the 
initial control \\p{0) — C(0)||x < S implies that for any t >0, \\p(t) — ^{t)\\x < rj. 
Otherwise, the solution ^ is said to he unstable with respect to the X norm. 

Of course, instability will then be interpreted as bad confinement, and stability as good 
confinement. 

We will investigate stability and instability around the following steady states, modeling 
the good and bad curvature sides: 

pf'-^x^) = (^''«'^'+ = 1 - ^, /i'"^'^'- = ^) (1.8) 

p9"°'^{xi) = (^p300d,+ ^ ^^^900d,- = 1 - ^) (1.9) 

which actually model a linear transition between the hot and the cold plasma. Indeed, for 
^bad^ the temperature of the plasma is given by: 

r(i,x) = T-^ + r+ (i-^) 

whereas for it is given by: 

T{t,oc) = T+^ + T- (l-^) 

Hence the profile of the temperature is a straight line with a slope equal to the so-called 
temperature gradient '^^^"^ . We mention that such linear profiles seem physically rele- 
vant in the edge of the tokamak (according to the graphs in [10] or [27]). 

Despite this rather rough model, our predictions will qualitatively correspond to obser- 
vations made by physicists. 

Our first objective will be to confirm the linear scenario exposed in the heuristic study 
by exhibiting a growing mode with maximal growth for the linearized operator in the bad- 
curvature area (but only when the temperature gradient "^^^"^ is not too large) and by 
showing that there is no such growing mode in the good- curvature area. 

Then our aim is to show that nonlinear instability also holds. As system (II. 6p looks 
a lot like 2D Euler, it is not so surprising that techniques allowing to pass from spectral 
instability to nonlinear instability for 2D Euler may apply here. On the topic of stability 
and instability of ideal plane flows, we mention some recent developments; let us neverthe- 
less emphasize that this list is by no means exhaustive. In |11) . Grenier proved instability 
in the velocity norm for some shear flows with zero Lyapunov exponent. In [3], Bar- 
dos, Guo and Strauss, following a method introduced in |12j , proved instability in the 
vorticity norm around some steady states. It is assumed that the linearized operator has 
a growth exceeding the Lyapunov exponent of the steady flow. We also mention the paper 
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of Vishik and Friedlander |25j where instability in the velocity norm is proved under 
the same assumptions on the steady states. The best result available by now is due to 
Lin |18) . Under rather general assumptions on the steady states (in particular, there is no 
assumption on the growth of the linearized operator), he showed nonlinear instability in 
the vorticity norm and in the same time, that velocity grows exponentially in the 
norm. In this work, we will obtain similar results to those of Lin. 

On the other hand, let us emphasize that in our linear analysis, when the tempera- 
ture gradient '^^^ exceeds a threshold, there is no growing mode in the bad curvature 
side. We will show that the nonlinear equations inherit this linear property. This rather 
unexpected stability phenomenon can be interpreted as the so-called High Confinement 
mode (H-mode for short), by opposition to the "standard" regime, referred to as the Low 
Confinement mode (L-mode for short). The H-mode is a high-confinement regime obtained 
by heating of the plasma and triggered when the heating power exceeds some threshold. 
It has been experimentally observed by physicists for a long time: it was discovered in 
the ASDEX tokamak [2^, we also refer to |10], [E] and ([27], Section 4.13). We can also 
remark that the H-mode is accompanied with an increase of the gradient of temperature 
(pJ]J, [27J). These experimental observations fit very well with our qualitative results. 

There exists a huge literature in physics on this particular topic. Nevertheless, despite a 
huge amount of works, the H-mode is still rather mysterious. Its understanding, especially 
the mechanism of transition from L-mode to H-mode is crucial for fusion research. To the 
very best of our knowledge, the H-mode has never been rigorously justified at the nonlinear 
or even at the linear level, with such a simple model. 

It is sometimes believed that the formation of confining transport barriers is due to a 
sheared E x B flow. In some sense, our model shares similarities with linear shear flows 
(the linearized equations are similar); thus our study is not in contradiction with these 
considerations. 

We finally mention that in the physics papers, the H-mode is most of the time numeri- 
cally investigated with more complicated models (including more physics, such as the effect 
of collisions, friction, energy sources), we refer to [6j and references therein. Our model 
can be seen as a two-temperature caricature of the model of [6]. The transition to the 
H-mode is also numerically investigated in |16) , where the existence of thresholds is shown. 
Unfortunately, we were not able to find any analytical formulae for those thresholds that 
we could have compared with ours. 

The main results proved in this paper (CoroUarv 15.11 and Theorem 16. ip are gathered in 
the following theorem: 

Theorem 1.1. For system il.b]) : 

1. (Nonlinear stability) 

The equilibrium ^9ood nonlinearly stable with respect to the norm. 

If the temperature gradient "^^T"^ satisfies: 



T+ -T 



1 




L 



then the equilibrium fj,^"''^ is nonlinearly stable with respect to the 

(Nonlinear instability) 

If the temperature gradient satisfies 



norm. 



T+ -T- 4 




L 57r2 
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There exist constants So,r]Q > such that for any < 5 < Sq and any s > there 
exists a solution p to U.6\) with ||p(0) — //''"'^||_f/s < 6 but such that: 

\\E(ts)\\L2>m (1.12) 

denoting E(ts) = VA~^{p'^{ts) + p'^its) — 1) the electric field at time ts = 0{\ log 6\). 
In particular, the equilibrium pf'"''^ is nonlinearly unstable with respect to the norm. 

1.3 Organization of the paper 

The present paper is organized as follows: section [2] is devoted to the study of the limit 
e — 7- for the system (|1.4p . In section [3] we present the simplified bi-fluid model we study 
in order to investigate stability and instability for the plasma. Section U] is dedicated to 
the study of the linearized system around the steady states p9°°'^ and particular 
we show the existence of dominant growing mode in the "bad curvature" region, provided 
that the gradient of temperatures is not too large. If the temperature gradient exceeds 
some threshold, then there is linear stability. In section O we are concerned with the 
nonlinear stability property for the "good curvature" region and for the "bad curvature" 
region for large enough temperature gradients (referred to as the high confinement mode 
in plasma physics), which will be achieved by exhibiting remarkable energies around the 
steady states. In section [6l for small enough temperature gradients we pass from linear 
spectral instability to nonlinear instability in the vorticity norm, using a high order 
approximation method introduced by Grenier. Then we prove that the electric field also 
grows exponentially in the norm, by using the energy exhibited in the previous section. 



2 Gyrokinetic derivation of the equations 

Fohowing Frenod and Sonnendriicker ([7]), we can use two-scale convergence tools in order 
to derive the gyrokinetic equation we are interested in. We shall not dwell on the rigorous 
derivation of this system since the justifications in two dimensions are essentially done in 

m- 

First of all, let us recall precisely the two-scale convergence notions (due to Nguetseng 
|21) and Allaire [T]) we will use in this section. 

Definition. Let X be a separable Banach space, X' be its topological dual space and (., .) 
the duality bracket between X' and X. For all a > 0, denote by Ca(M, X) (respectively 
La (M; X') ) the space of a-periodic continuous (respectively L'^ ) functions on M with values 
in X. Let q S [1; oo[. 

Given a sequence {ue) of functions belonging to the space L'^ (0, t; X') and a function 
U°{t,9) £ L'^'{0,T-Li{R-X')) we say that 

Ue 2-scale converges to 

if for any function ^ G ^"^(0, T; Cq(M, X)) we have: 

lim^ (^u,{t),'i' (^t,-^ dt^ = ^ {U'^{t,T),^{t,T)) drdt (2.1) 

The new variable r has to be understood as a "fast-time variable" which describes 
the fast time oscillations. As for weak-star convergence in L^ spaces, one can show that 
boundedness implies 2— scale convergence in spaces. 



7 



Theorem 2.1. ffHf. l^) 

Given a sequence (m^) bounded in L'^ (0, there exists for all a > a function 

C/° G L«'(0,r;L^(M;X') such that up to a subsequence, 

2-scale converges to 

The profile is called the a-periodic two scale limit of and the link between and 
the weak- * limit u of is given by: 

-ruldT = u (2.2) 



For the reader's sake we recall the main arguments and refer to [T] for the complete 
proof. 

Sketch of proof. Let a > 0. We can consider (p^^ '■ 

cT 



^ € L'^{0,T;CaiR,X)) ^ u,{t)^ (^t, 



dt 

e 



and show that it is a continuous linear form on L''(0, T; Cq,(M, X)), so that it can be 
identified with a unique 11^ in L''(0, T; Cq,(M, X))'. Then we can show that 17^ is uni- 
formly bounded in ^9(0, T; Ca(M, X))' ; thus, since L'?(0, T; C„(M, X)) is a separable Ba- 
nach space, weakly-* converges up to a subsequence to some U in L'?(0, T; Ca(M, X))'. 
Using Riesz's representation theorem, one can show that it can be identified with some 
G L'^ {0,T;La (M;X')) and that two-scale converges to Uq. 

□ 

We can now state the main result of this section: 

Proposition 2.1. For each e, let f^ be a global weak solution to ljil-4\ ) in the sense of 
Arsenev . 

Then, up to an extraction, f^ 2-scale converges to a function F: 

F{t, T, X, v) = G{t, X + n{T)v, R{t)v) (2.3) 

and G satisfies: 

8.0 + (a nrm, r, . + 7J(-r)„)<ir + (^,^(^_ I^tl"! + „|,) ) ^ V.G 

+ [h Jo' R{r)£{t,r,x + W-r)v)dr + (_"jj^j'^ +/^^))) -V.O = 

G\t=o = fo 
E = -V^V 

-Ay = j G{t,x + n{T)v,R{T)v)dv -1 
denoting by R and TZ the linear operators defined by : 



(2.4) 



Rir) 



cos T — sm T 
sin r cos r 



, 7^(T) = (i?(-7r/2) - i?(-7r/2 + r)) . 
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Proof. We do not wish to develop the very beginning of the proof since it is strictly identical 
to the one given in [?]. 

The first step consists in deriving the so-called constraint equation. To this end, let 
^{t,T,x,v) be a 27r-periodic oscillating test function in r and define: 



Then we can write the weak formulation of the Vlasov equation against and pass to 
the two-scale limit. We find that the two-scale limit of f^{t,x,v), denoted by F{t,T,x,v), 
satisfies the following equation: 

drF + v_L.VxF + v Aez.VvFa = 0, (2.5) 

As a consequence, F is constant along the characteristics meaning that there exists a profile 
G with: 

F{t, T, X, v) = G{t, X + n{T)v, R{t)v) (2.6) 

where R and TZ are defined in the proposition. 

The next step is to determine the profile G. We introduce the filtered function g^: 

g,{t, X, v) = Mt, X + n{-t/e)v, R{-t/e)v) (2.7) 

which represents the number density from which we have removed the essential oscillations. 
Notice that this function is chosen so that two-scale converges, as well as weakly-* 
converges to G. 

We easily obtain the equation satisfied by g^^: 

dtge+n{t/e)E,{t, X + n{-t/e)v).V.xg, 
+R{t/e)E,{t, X + n{-t/e)v).Vyge 

-n{t/e) ((x + n{-t/e)v)^.{R{-t/e)v)^) .V^g, (^.S) 
-R{t/e) ((x + ni-t/e)v)i.iR{-t/e)v)^^ .V,g, = 

We now pass to the limit in the sense of distributions. We can prove that the following 
convergence holds for the nonlinear terms (using elliptic regularity for the electric field to 
gain some compactness): 

n{t/e)E,{t, X + lZ{-t/e)v).Vxg, [ n{T)£{t, r, x + lZ{-T)v)dT.V ^G (2.9) 

1 /■^'^ 

R{t/e)E,{t, X + n{-t/e)v).V^ge — / R{T)£{t, t,x + n{-T)v)dT.V^G (2.10) 

2vr Jo 



Likewise, we have the following convergences for the last two terms (here there is 
basically nothing to justify since these are linear quantities): 

-n{t/e) (^{x + n{-t/e)v)i.{R{-t/e)v)^^ .V^g, 

j^J n{T){{x + n{-T)v)UR{-r)v)^)dT.VxG ^ ^ 



9 



-R{t/e) [{x + n{-t/e)v)i.{R{-t/e)v)^) .V,g, 

1 i-2n (2.12) 

— / R{T){{x + n{-T)v)i.{R{-T)v)^)dT.V,G 

JO 

We then compute the foUowing quantities: 

1 Z"^'' N// , -n/ N ^ /r./ N ^±^J V2(-Xi+V2) 



R{r){{x + TZ{-r)v)i x (i?(-r).)^)dr = (2-13) 

- ^"7^(r)((x + 7^(-T)^;)l X (i?(-T)t;)^)dr = f , + ^2) \ ^3.14) 

i^io ^ V ;n V V \^^2(a;i-^;2)-i(^^f + ?^i)y 



27r 

27r 

This concludes the proof. 



□ 



Qualitative interpretation of the gyrokinetic system 

The influence of the variations of B is given by the drift/acceleration terms: 

-Vi{v2 -Xi) \ ^ r ^( ^2(-2;i +V2)\ ^ 

V2{XI - V2) -\{vl+ VD) -^^^ + + ^2) ) 

Let us imagine that there is no electric field in the asymptotic equation (12.40 . Then, 
the characteristics are given by the following ODEs: 



dx _ ( -Vl{v2 - Xi) 



\V2{XI - V2) - U^l + V^) 



dv 
dt 



V2{-Xl + V2) 
-Vi(-Xi + V2] 



(2.15) 



At first sight, this dynamical system seems a bit complicated with some unpleasant 
quadratic terms. Actually, this system has some nice invariants. 
First, notice that 

— (xi -V2) = 0. 

This means that xi = V2 + Ci (with Ci = xi{0) — ^2(0)). The equation for v can now be 
written in the simple form: 

? = (2-16) 
dt V C^i^i J 

The velocity is thus periodic (and we could compute it easily). Notice also that 

|(.? + .i) = o, 

so that vj + vl = C2 (with C2 = Vi(0) + ^1(0)). 

We get as well a periodic motion for xi (since xi = V2 + Ci) . Finally we notice for X2'- 

j^{x2 + vr) = -\{vl + vl) = -\c2 (2.17) 

The motion along the 62 direction is hence a sum of a periodic motion plus a fall which 
only depends on the initial velocity of the particles (and not on their position). Such a 



10 



100 



200 



300 



400 



500 



600 



700 



800 



900 



1000 



Figure 2: X/B drift (in the x-axis: 
time and in the y-axis: X2) 



Figure 3: Motion of a particle starting at (0,0) 
in the shce of the tokamak (in the x-axis: xi 
and in the y-axis: X2) 



drift of the particles "in the bottom" of the tokamak and depending only on the square 
of their velocity is predicted by physicists and is often referred to as the VB drift ([ST], 
Section 2.6). To support our discussion we give some graphs of the characteristic curves 
(figures [2] and [3| . 

Likewise, Frenod and Sonnendriicker introduced in \J\ the new variables: 



the so-called Larmor radius variable. With these, they showed that the terms in (12. 4p 
involving the electric field were qualitatively close to the drift 



which corresponds to the gyroaveraged electric drift (|2?]. Section 2.11). 

3 The simplified mathematical model 

To be completely rigorous, system (12. 4p is the system we would have to study in order 
to investigate good or bad confinement. Nevertheless, at least at first sight, its algebraic 
structure seems to be too complicated. Morevore, it is not so clear how to choose a steady 
state describing the physical situation we want to study. 

Consequently we will make several approximations (some of them being quite rough) 
on (12. 4p in order to get a more tractable model. 

3.1 A drift-kinetic system 

A first step is to obtain a simplified kinetic system, whose dynamics is close to system (j2.4p . 
We therefore consider the following drift-kinetic equation, which is actually a classical 




the so-called guiding center variable and 



w = —V 
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physical model ([2?J, Section 2.11). It is commonly used for numerical simulations (see for 
instance the GYSELA code M): 



dtf-^\v\^d^J + E^.V^f = 
-A,V = I fdv-l 



This system can be heuristically derived from Newton equations with some elemen- 
tary physical considerations (see for instance |27] . Section 2.6); unfortunately we were not 
able to derive it rigorously from (II. 4p or (12.40 . Nevertheless, considering the qualitative 
study of last paragraph, this seems to be a reasonable model, if we make the following 
approximations : 

• We neglect the oscillations in time, which amounts to get rid of the explicit depen- 
dance on the fast time variable r for the electric field. This can be justified if we 
consider well-prepared initial data: we refer to the work of Bostan [1]. 



We neglect the gyroaverage operators: 
27r 



1 f'^^ 

— I E-^(t,Xc + R(T)w)dT E^it, Xc 
!7r Jo 



which is reasonable if we consider that the variation of the electric field across a 
Larmor radius is negligeable. 

3.2 The bi-temperature drift-fluid system 

In order to get a simplified fluid model, we assume that the plasma is made of two phases, 
one being the cold plasma (with low velocities, low temperature T~ and density p~) and 
the other the hot plasma (with large velocities, large temperature and density /9+). Of 
course, we take < . 

Hence, we assume that the solution to (13.11) takes the form: 



f{t,x,v) = p+{t,x)u+{v)+p-{t,x)v'{v) (3.2) 
where p"*" (resp. p") is a positive density such that the total mass is equal to 1, that is: 



/ 



(p^ + p )dx = 1. 
Furthermore, and are measures defined by: 



1 

"W2f^^l^l=^ 

Considering that transverse particle velocities are isotropically distributed is physically 
relevant for such magnetized plasmas, as indicated in [23j. 

We observe that J dv^ = 1 and J vdv^ = 0. Thus, the charge and current densities 
are given by: 

p(t, x) := J fdv = p'^{t, x) + p~{t, x) 
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and 



u{t, x) :- 



0. 



In addition, we have: 



/ fvdv 
JJdb 

v\^dv^{w) 
vf'dv^ {w) 



1 
2 

2 

The kinetic temperature T(t, x) of the plasma is then given by: 

2 p{t,x) 
_ p+{t,x)T+ + p-{t,x)T- 
p+{t,x) + p-{t,x) 



(3.3) 



We moreover assume we can decouple the transport equations satisfied by p^ and p . 
We get in the end the macroscopic system: 

( dtp+ - T+d^,p+ + E^.V^p+ = 
dtp- - T-d,,p+ + E^.V^p- = 

E = -V^V (3.4) 
-A^V = p+ + p- -I 
. ip^,P~)\t=o = {po,Po) with /p+ + pg = 1 

with x G [0,L] X R/LZ. 

As noticed in the introduction, this systems looks like 2D incompressible Euler, but 
with two kinds of vorticities. 

Here, the constant L > stands for the size of the box . The periodicity with respect 
to X2 is physically justified if we consider that L is small enough with respect to the size of 
the tokamak, so that we can decompose it in many identical cells of size L (see Figure 1). 

We now have to impose some relevant boundary conditions on xi = 0, L: 

• For the Poisson equation, we opt for the perfect conductor assumption on xi = 0,L 
(which is the ideal case for plasma physics models). 



(3.5) 
To this end, we can impose the following Dirichlet boundary condition on xi = 0, L: 

V = (3.6) 



From the fluid mechanics point of view, we observe this corresponds to the classical 
no slip condition. 

• For the transport equation, we actually do not need any boundary condition. There 
is indeed no entering or leaving trajectories, since the linear "drift" operator only 
entails a motion along the 62 direction, and = on the boundaries xi = 0, L. 

Following classical works on the Cauchy problem for the 2D incompressible Euler sys- 
tem (we refer for instance to the book of Majda and Bertozzi |20)). we get the following 
global existence and uniqueness result of strong and weak solutions to (j3.4p : 
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Theorem 3.1. Let po = (Po ,Po ) ^ (^H(0>-^^) x ^/LZ)f with Pq , Pq non-negative and 
HPo +Po)dx = 1. 

1. (Kato, /iTj/) If Pq is H'^ (with s > 1) then there exists a unique classical solution p 
to [T^ inCf{[0,oo[,H'')nCl{[0,oo[,H''-^) with initial data pQ. 

2. (Yudovic, f28^) If po G then there exists a unique global non-negative weak 
solution p € Lf{L^ n L°°) to {3.4^ with initial data pQ. 

Sketch of proof . 1. The existence of a global strong solution follows from a fixed point 
argument, as in the classical work on 2D Euler by Kato [17]. Actually, dealing 
with two "vorticities" only slightly modifies the main lines; therefore, for the sake of 
brevity, we only recall here the main arguments of the proof. 

We may consider the map F : {(,^,^~) i— )■ {p^ , p^) where {p^,p^) is solution to: 

dtp+ + {E^ - T+ea) .V,p+ = 
dtp- + {E^ -T-e2).V^p- =0 

y = on xi = 0,L 

{p'^,p-)\t=o = {po,Po) with. Jp+ + p- = l 

which is well-defined if is smooth enough thanks to the characteristics' 

method. As in Kato's proof, for any T > 0, one can show that F is continuous 
on some convex compact S of (C([0,r] x [0,L] x R/LZ) f and that F{S) C S. The 
existence of a fixed point is finally a consequence of Schauder's theorem. The crucial 
points are: 

• Establishing some log-lipschitz estimate on the electric field (in this case with 
a constant involving the norms of the two vorticities), which is obtained 
exactly in the same way as for Euler's equation. 
If £ L°° and E is solution to the elliptic problem 

E = -V^V 

-A,y = e+ + r-i (3.8) 

y = on j;i = 0, L 

then there exists C depending on L but independent of such that for all 

(x,y) G ([0,L] X M/LZ), we have: 

\E{x) - E{v)\ < C {Wth^ + liriU-) \x - y\ log+(|x - y\) (3.9) 



where log^(s) 



1 — log s if s < 1 
if s > 1 



• The L°° norm of each "vorticity" p+ and p is conserved by the transport 
equations, since div^; (^E-^ — T^e2) = 0. 

Finally uniqueness is obtained by an energy argument, exactly as in Kato's proof. 

For the bi-dimensional Euler equations in a general domain, the result is due to 
Yudovic |28) . His main arguments can be easily adapted and reproduced in our 
case, with a crucial use of the two above points. The result also follows by a simple 
adaptation of the alternative proof given in (|20). Theorem 8.1), which consists in 
regularizing the initial data, solve the smoothed Cauchy problem and then pass to 
the weak limit thanks to a crucial compactness result (pQ], Proposition 8.2). 

□ 
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3.3 Modeling of the plasma equilibria 

We now consider the following steady states, in order to model what is happening in the 
"good curvature" or the "bad curvature" side of the plasma, the only difference being the 
relative position between the hot and the cold plasma: 

• in the "bad curvature" region: 



hot plasma 



L 



cold plasma / 



(3.10) 



in the "good curvature" region: 



good 



(Xl) 



good,- 



^1 i.good,- 



L 



(3.11) 



These are steady states of (|3.4p and the associated electric field is zero. 
We observe here for ^9°°'^ temperature T{t, x) is: 



T{t,x) 



L 



1 



L 



(3.12) 



Such linear transitions between the cold and hot plasma are the most simple model one 
can think of. We observe that the slope of the line is equal to "^^^ which is referred to 
as the temperature gradient in this paper. 



We now investigate stability and instability for fj'"''^ and 11^°°'^. 



4 Linear instability in the "bad curvature" region 

We first consider the case of the "bad curvature" region, for which we expect to obtain 
instability. The equilibrium writes: 



1 _ £1 
^ L 

L 



The first step before trying to prove any instability property for the nonlinear transport 
equations consists in investigating the problem of instability for the linearized operator 
around /i''"'^. We accordingly consider the following linearized system (for t > 0, x G 
[0,L] X M/LZ): 



( dtp+ - T+d,,p+ -^ = 



dtp- - T-d,,p- + ^ = 

E = -V^V 

-A^V = p+ + p- , y = on xi = 0, L 
, (p^,P~)|t=o = (Po ,Po ) witli /Po +Po =0 



(4.1) 



4.1 Looking for unstable eigenfunctions 

We look for special solutions under the form /9fc(t, x) 
can be restated as: 
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f dtih+{t)gk{x)) - T+d,,iht{t)gk{x)) - f = 
dtiKit)9kix)) - T-d,,{h-{t)gk{x)) + f = 

E = -v^y 

-A^y = {hl{t) + hl{t))gk{x) , y = on XI = 0,L 



(4.2) 



We take g^ with the particular form gk{x) = sin (-^7rxi)e* ^^T^^ (^-^^^^1^ ki,k2 G Z*), so 
that is an eigenfunction for the laplacian with the considered boundary conditions. It 
satisfies indeed: 

1,2 ■ 



^9k{x) 



-TT 



L2 ^ L2 



(4.3) 



and also ^(fc = on xi = 0, L. 

The solution to the Poisson equation is then given by: 



and thus we have 



ihkit) + h+{t))gk{x) 



-i2k2/L 



-{Kit) + ht{t))gk{xy, 



7r(§+4§) 

which leads us to study the following first order ordinary differential equation: 

\ 

(htit)\ 



+ 1/L 



i2k-2 



-(4+4) 



-2mT-k2 - 



»2fc2 



-(4+4) / 



(4.4) 

We want to compute the eigenvalues of the matrix; its characteristic polynomial states: 



X"^ + 2mk2 {T+ + T-)X- ATr^k^T+T' - 



Akl 



and its discriminant: 



A = -47r^A;i(T+-T-)2 + 



IGk^L 
kf + Aki 



-4Ai(T--T-)((T--T-)--,^^i^) (4.5) 



We recall that by definition, 

T+ -T~ >0 
We can now distinguish between two cases: 
• First case: 



4 T+-T- , , 

> —[— («> 

In the case where the gradient of temperature is not too large, then there exist fci, /c2 
such that A > 0. We consequently obtain two complex roots, one of which has a 
stricly negative real part equal to — In other words, this shows the existence of 
an unstable mode. 
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• Second case: 



4 T+ 
5^ - T 



(4.7) 



In the opposite case, we always have A < and consequently, we are not able to find 
a growing mode ! 

This phenomenon may at first sight look like a mathematical artifact due to the peri- 
odicity constraint in the X2 direction. Nevertheless as explained in the introduction, 
the existence of such a threshold is well known in plasma physics, beyond which one 
can expect tremendous confinement properties. In very rough terms: heating brings 
stability. The stable mode is referred to as the H-mode, by opposition to the L-mode. 

Remark 4.1. In the "good curvature" region, that is around figood: the linearized system 
states: 



dtp 
dtp 
E = 



T-d,,p- 



L 

Ml 

L 








p^ + p~ , y = on xi = 0, L 
With the same method, we obtain the following ordinary differential equation: 



(4. 



dt 



Kit) 
Kit) 



/-2mT+k2 



i2k2 



+ 1/L 



L ^ L 



i2k2/L 



i2k2 



L ^ L 



-2mT- +k2- 



We consequently have to look for the roots to the polynomial: 



i2k2 



Kit) 
Kit) 



(4.9) 



X'^ + 2mk2v'^X - An'^T+T-kj + 



Aki 



fr2 i.2 



(r+ -T-) 



In this case, one can check as before that the discriminant is always stricly negative, so 
that the roots always have a vanishing real part. As a result, we do not find any unstable 
mode by this method. Note that we only consider this fact as a good and encouraging 
indication for stability around p^°°'^. Actually, we will never use it when proving nonlinear 
stability in section [5j 

Remark 4.2. We finally mention that the quest for unstable modes seems more difficult 
in the kinetic case, since one has to deal with the continuous velocity space. A very famous 
criterion in the Vlasov-Poisson case was given by Penrose [22] and rigorously studied later 
on by Guo and Strauss [13] . 



4.2 On the spectrum of the hnearized operator around /i''"'^ on H'^([0, L] x 
M/Z), s > 

We assume here the existence of unstable modes for the linearized operator around p^"''^, 
that is in the situation where we have: 



4 T+ -T- 
5^ ^ I 
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The main tool we use now is a variant of a classical theorem by Weyl, stated for 
instance in the paper of Guo and Strauss |13) and proved by Vidav in |24) . Basically, 
it gives informations on the spectrum of some compact perturbation of a linear operator 
which has no spectrum in the half-plane {Rez > 0}. It entails the existence of a dominant 
unstable eigenvalue provided the existence of at least one growing mode. 

Theorem 4.1 (Weyl). Let Y be a Banach space and A be a linear operator that generates 
a strongly continuous semigroup on Y such that \\e~^^\\ < M for all t > 0. Let K be a 
compact operator from Y to Y . Then (A + K) generates a strongly continuous semigroup 
^-t{A+K) and a{—A — K) consists of a finite number of eigenvalues of finite multiplicity in 
{Re A > 5} for every 5 > 0. These eigenvalues can be labeled by: 

ReAi > ReA2 > ...ReAiv > S 

Furthermore, for any 7 > ReAi, there exists some constant C-y such that 

||e-*(^+^)||y^y <C^e*^ (4.10) 

Corollary 4.1. Let s > and 

Y = {yi,y2e H'{[0,L]xR/Zf, j{yi+y2)dx = 0}. 
Let M be the linear operator defined by: 

with E2 = —8x2^) —Ay = {g^ + g~) and with V = {) on xi = 0, L. 

Then there exists an eigenvalue A with a non-vanishing and maximal real part associated 
to a C°° eigenvector. Furthermore for any 7 > Re A, there is a constant 0(^,3) such that 
for all t>0: 

\\e-'^'\\Hs^Hs<C{j,s)e'^ (4.12) 



Proof. The linear operator A, defined by 

A:g^ 



-T+dx,g^ 
-T'dx^g' 



is clearly an isometry on Y (indeed we know how to explicitly solve the semi-group). The 
operator K is defined by 

/-Ml' 
K-.g^i^^ 

This operator is compact on Y thanks to standard elliptic estimates. 

Moreover, we have shown in the last paragraph the existence of an unstable eigenfunc- 
tion for the linearized operator that belongs to any i/'^. We can therefore apply Weyl's 
theorem which gives the existence of an eigenfunction associated to an eigenvalue with a 
non- vanishing and maximal real part. At last, the estimate in the corollary follows directly 
from the estimate given in Weyl's theorem. 

□ 
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In the following, we denote for any h & with J hdx = 0, A the unique solution 
u in to the problem: 

-An = h 
M = on xi = 0, L 

We give now a lemma which tells us that any eigenvector associated to a non vanishing 
eigenvalue for the linearized operator any Lp' actually has regularity. 

Lemma 4.1. Let p = (p"*", p~) G (I/^)^ with J (p"*" + p~)dx = and A / such that: 

-r+a,,/9+-^a,,A-i(p+ + p-) = \p+ 
-T-d^,p- -^d^,A-\p+ + p-) = Xp- 

then {p-^,p-) G C°°([0,L] x M/Z) 

Proof. The principle of the proof is to show by recursion that p G H'^, for any G N*. 

For = 1, we can observe, thanks to elliptic estimates, that dx2P G L^. Indeed, we 
have the identity: 

- T+d,,p+ = ^a,,A-i(p+ + p-) + Xp+ (4.13) 

Hence, dx^p^ G L^. Likewise, dx2p~ G L^. 

We can apply the differential operator dx^ to the equation satisfied by which entails: 

-T+dx,dx2P+-^dx,dx2A-\p+ + p-) = xdx,p+ 

Then we multiply by dx^p^ and integrate with respect to x: 

^ll^a;iP"^|li2 = J -T^dx^dx^P^dx^p'^dx - J J dxj^dx^A'^ip'^ + P~)dxiP'^dx 
Thanks to the periodicity with respect to X2, we get: 

J dx^dx2P^ dxiP^ dx = 1/2 j dx^ldxiP'^^ dx = 
Then using Cauchy-Schwarz inequality: 

M\dxiP'^\\l2 < ^Wdx.dx^A'^ip'^ + p~)\\L^\\dxiP'^\\L2 
As a result we showed that: 

X''\\dx,p+\\L^ < ^\\dx,dx2A-\p+ + p-)\\i^2 (4.14) 

By standard elliptic estimates the right-hand side is finite since p'^ et p~ belong to L^. As 
a result we have proved p e H^. 

We can then conclude by recursion. Let us assume that p G H'^, for some A; G N*; we 
prove that p G H^^^. 

Let a,/3 G N such that a + (3 = k. We set = d^^dx^- Then d^p belongs to and 
satisfies the equation: 

-T+dx2dkP+ - ^dkdx2A-\p+ + p-) = XdkP+ 
-T~dx2dkp~ - 2^kdx2^''^{p'^ + P~) = ^dkp~ 
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By elliptic regularity, 9^5x2 A ^{p^ + p ) G H^. Thus, we are in the same case as for 
L'^ ^ H^, which entails that p € H^^^ . 

□ 



5 On nonlinear stability 

Let p = ^ solution to the nonlinear transport equation p.4p . that we recall here: 

' dtp+ - T+d,,,p+ + E^.V^p+ = 

dtp- - T-d^,p+ + E^.V^p- = 

E = -V^V 

-A^V = p+ + p- -1 

y = on xi = 0,L 
, (/5+,P")|t=o = (Po ,Po ) with + =1 

We begin with a very simple observation in the limit case T = = . In this 
situation, setting p = p'^ + p~ , p satisfies the usual 2D Euler equation in vorticity form, 
with some linear drift term: 

dtp + E^.V^p-Td^,p = 
=V^A-\p-l) 

We investigate stability around the steady state p, = p^ + p~ = 1. It is well-known that 
any norm of the "vorticity" p — 1 is non-increasing, that is: 

\\p{t) - < \\m - IIIlp (5.1) 

which clearly entails nonlinear stability. 

The idea in the general case T+ > in order to show nonlinear stability is to obtain 
some similar nice energy estimate. 

Theorem 5.1. For any inital data pQ G with J{Pq + P())dx = 1, the solution p to 
satisfies the following statements. 

• Around the "good-curvature" steady state, the following functional is non-increasing: 

£{t) = \\p- p^^^Xl^ + j I vypdx < £:(0) (5.2) 

with \\p — /i||^2 = ll/o^ — ^^'''11^2 + \\p^ — ^^11^2 and W = VA^^(p+ + p^ — !)• 

• Around the "had- curvature" steady state, the following functional is non-increasing: 

m = \\p- i^'^'wh - ^T+^_T~) j i^^i'^^ ^ -^(0) (^-^^ 

Remark 5.1. We observe that in the functionals, the first term corresponds to the enstro- 
phy in fluid mechanics (in a modulated form adapted to our needs). The second term can 
be interpreted as the kinetic energy of the fluid (whereas from the plasma physics point of 
view, this is the electric energy). 
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As an immediate consequence of this theorem, we obtain stability, in the "good- 
curvature" side, but also in the "bad-curvature" side, for large enough temperature gradi- 
ents, like for the linearized equations. 

Corollary 5.1. The equilibrium ^^ood nonlinearly stable with respect to the I? norm. 
If the temperature gradient satisfies 

—r^ > -2 (5-4) 

then the equilibrium ^'""^ is nonlinearly stable with respect to the norm. 
Proof of the corollary. Thanks to the energy identity (15. 2p and to the Poisson equation: 

Wp-^'^'^'Wh < \\P - P'°"'\\h + ^(^+^_ / \^V\'dx 

< MO) - ^^-ii. + I ivvm'dx 

< MO) - /.^-ii. + C ^^^^\^_^ MO) - M^-'^lli. (5.5) 

This means that for any 7] > there exists 6 > such that if ||p(0) — H^°°'^\\l2 < 6 then for 
any t > 0, \\p — ij,^°°'^\\]^2 < r]. In other words, ^9°°'^ is nonlinearly stable for the norm. 

Around the bad-curvature steady state, we have shown that the following quantity is 
non-increasing: 

Ht) = \\P - p'^'Wh - rfrr^^ / l^^l'^^ (5-6) 



L{T+ - T 

We can easily prove with the help of Fourier variables the Poincare-like inequality: 



\VV\'^dx<^\\p- ^J^'"^\\l^ (5.7) 



So we get: 



Hence, 



1 - \\p - p^-'Wl^ < ||/>(0) - p^'^Xl^ (5.8) 

As a consequence, there is Lp' nonlinear stability in the bad curvature side, provided 
that 

T+ -T- 1 



> 



L 7r2' 
Otherwise, we can not deduce anything. 

□ 

Remark 5.2. Note also that in the linear discussion, there was stability provided that 

T+ -T~ 4 

^7~T— |- rri — r A "In 

We do not know what happens for — j — G [5;;^' '^l ^^^^ nonlinear case. Maybe one 
could expect to observe some bifurcation phenomenon. 
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Remark 5.3. Actually the decrease of £{t) tells us a little more than just stability. 
Indeed, there exists C > such that for any 5 > 0, if ||/9(0) — /^^°'"^||^2 < then for any 
t > 0, 



This means in particular that for large values of — T , better confinement is obtained, 
which is qualitatively in agreement with experimental observations (|10)). 

We now prove Theorem 15.11 We first give a technical lemma in which will help for the 
proof. 

Lemma 5.1. For /j, = ^sood ^bad^ have for any t > 0, 

j E2 (p+ -fj.+)dx = - j E2 {p- - fi-) dx. (5.10) 
Proof. In order to prove this identity, one can simply compute: 

J E2 {{p- - fi-) - ip+ - dx = j E2{p'^ + p- - iJi^ - p- -2{p+ - iji+))dx 

= j E2 {AV - 2{p+ - p-^)) dx 
= -2 j E2 (p+ - dx 
Indeed, thanks to periodicity with respect to X2 and since 8x2^ = on xi = 0, L, we get: 
J dx^VAVdx = - j dx.,VV.VVdx + j diw{dx.,VVV)dx 



vv 



dro — - — ] dx = 



2 

If we make the same computation, by symmetry, we can also observe that: 
J E2 {{p- - p-) - (p+ - p+)) dx = 2 j E2 {p- - p-) dx 



□ 



Proof of Theorem \5.1\ We will only focus on the proof of the conservation of £{t), the 
proof being very similar for J-{t). For the sake of readibility, we write p instead of p3"°'^ 
until the end of the proof. 

We observe that the equation satisfied by [p — p) reads in this case: 

and E = —VxV with —AV = p^ + p^ — p'^ — p^ and with Dirichlet conditions on the 
boundaries xi = and xi = 1. 

Taking the scalar product with (p — p) in the transport equation and integrating with 
respect to x entails: 

i\\P-t^\\h = l -T^P"- - P^)d- + / f (P" - (5-12) 
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Indeed, thanks to the periodicity with respect to X2^ we first have: 

j - - ix+)dx = j \d,,{p^ - ii+fdx = 

Similarly we have: 



(p - /" )(P - )dx = Q 
In the same fashion, with Green's Formula, and using div = 0, we have: 



since E2 = —8x2^ = on xi = 0, L. 
Now, by Lemma IS.ll we get: 



^ I dW{E^ip-pf)dx 

]-( f E^{p-fif.{-ei)dx2+ [ E^{p- pf.eidxA 

^ \Jxi=0 Jxi=l J 



^llp-^llia = J -^ip'^ - ^'''^)dx + J ^{p - p. )dx 

= -2j^{p^-p^)dx{=2j^{p~-p-)) 



We have: 



- j E2 (p+ -p^)dx = - j Vdx2 (p+ - dx + j div(y (p+ - p+) e2)dx 

V ' 

=0 

= ^J-V (^dtip^ - p+) + E^.Vxip^ - + §) dx 

= {-dtip-^ + p--p+- p-) + E^.Vxip-^ + p--p+- p-)) dx 

+ ^ J -T-Vdx.ip- - p-)dx 

= ^ Jv (dtAV - E^.Vx^v') + ^ / ^2(P" - p-)dx 

where we have plugged in the equation satisfied by {p~^ — by (p~ — p~), and also 
plugged in the Poisson equation. 
Finally we have: 

(1 - ^) - /^lli^ = / ^ {dtAV - E^.VxAv) dx 



So that: 
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To conclude the proof, we observe, using the Dirichlet boundary conditions and the 
periodicity: 

VdtAVdx = -~ ( [ \VV\'^dx] (5.13) 



dt2 



and 



/ 



VE^.VAVdx = j Vdiw{E^AV)dx 

= - j VV.E^AVdx + j A:w{VE^AV)dx 



The first term is equal to zero since E.E-^ = 0, the second one thanks to the boundary 
condition on xi = 0,L and to the periodicity with respect to X2- 
As a result we have proved that 

Actually the computations we have made are rigorously valid only for smooth solutions 
to (13. 4p . Nevertheless, these can be justified by smoothing the initial data, and then passing 
to the weak limit, which entails (|5.2p . □ 

Remark 5.4. We observe that we have also proved that: 
which yields: 

\\p''{t)-^^^\\l2 - \\p-{t) -^^-\\l2 = ||p+(o) - ||p-(o) -/i-lli^. (5.15) 

Remark 5.5. Let us add that the explicit form of the equilibria is crucial in the proof of 
the theorem. It would not work similarly if we had taken an equilibrium of the form: 

^^'''^=[l-^x,) 

with $ a smooth function. In this case it should be maybe more relevant to use the general 
Lyapunov functionals method of Arnold [2j. 

Likewise one can notice that the proof would have not worked if we had chosen any 
other boundary condition than Dirichlet. 



6 On nonlinear instability 

What we intend to show now is a property of nonlinear instability in the "bad curvature" 
region when the physical parameters satisfy: 

T+ -T- 4 

< 



L 57r2' 

This can be interpreted as a bad confinement property. We first recall that the equi- 
librium in this case is the following: 
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) 



Thanks to the existence of an eigenvalue with maximal positive real part for the lin- 
earized operator around /z''"'^, we can prove a nonlinear instability result. Indeed, using the 
method introduced by Grenier we are able to pass from the linear spectral instability 
to the nonlinear instability in the norm. Grenier's method was originally used to prove 
instability in the velocity norm for Euler ; we show here that this technique can also be 
adapted to show instability in the vorticity norm. 

The drawback of this method is that it requires high regularity on an eigenfunction 
associated to the dominant eigenvalue, which could be difficult to check in more complicated 
cases. In our case, we were able to prove such a smoothness in Lemma l4.ll 

For the sake of readability we will write ^ instead of ^^"''^ since there is no risk of 
confusion. 

Theorem 6.1. There exist constants 5o,r]i,rj2 > such that for any < 5 < 60 and any 
s >0 there exists a solution {p,E) to Ii3.4\ ) with \\p{0) — pWh" < ^ but such that: 



with ts = 0(1 log 

In particular, p is unstable with respect to the norm. 

Remark 6.1. This instability result is complementary to the stability result proved in 
Corollary 15. H since they involve the same norms. 

Remark 6.2. This instability result can also be obtained by techniques similar to those 
used by Bardos, Guo and Strauss in |3j for the 2D incompressible Euler system, which 
consist in proving that the solution p remains "close" in some norm to a growing mode 
associated to the maximal growth rate of the linearized operator M. 

Using some bootstrap argument introduced by Bardos, Guo and Strauss [3j, and a idea 
of Lin [18] consisting in: 

• estimating ||iiJ(t5)||j;^2 with Duhamel's formula 

• then studying the semi-group e~*^^ (where M is the linearized operator) on H~^, 

we can also prove the exponential growth of the norm of the electric field. 

We refer to the paper of Lin [19J which could be adapted to our case with some minor 
modifications. 

Here we will provide an alternative proof by using Grenier's method (which consists in 
proving that the solution p remains "close" in some norm to the growing mode plus some 
high order correction) to prove (16. ip and finally by using the energy of Theorem 15.11 to 
prove (|6.2p . 

Proof. We begin with some preliminaries on the linearized operator. Using the same no- 
tations as in paragraph 14.21 we consider M the linearized operator around p on 



\\p{ts)-p\\L^ > m 



(6.1) 



and: 



\\E{ts)\\L2>m 



(6.2) 



Y = {yi,y2 G H'{[0,L] x M/Zf, / {yi+y2)dx = 0}, fors > 0. 
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For s = 0, by Corollary 14. H we know the existence of an eigenfunction R associated 
to an eigenvalue A with maximal real part Re A. In addition, by Lemma |4.H R belongs to 
any H'^, s > 0. For the sake of simplicity we will assume that R is real and associated to 
the eigenvalue Re A. In the general case, since the linearized operator is real, the conjugate 
of A is also an eigenvalue so that one can consider by linearity real-valued growing modes 
and the following of the proof remains the same. 

We recall also that by Corollarv 14. H for any 7 > Re A, there is a constant C(7, s) such 
that for all s > 0: 

\\e-'^\\Hs^Hs<C{^,s)e'^ (6.3) 

Basically the idea of Grenier is to construct a high order approximation of the nonlinear 
equation, that is a more precise approximation than the "usual" linearized equation. Indeed 
instead of showing that / — /i is close to a well chosen eigenfunction, we show that it is 
close to the high order asymptotic expansion: 

N 

p\Sl=5u^+Y.5'u, (6.4) 

i=2 

where ui = ReP"'^^^ . Note than for any s > 0, we have: 

\\ui\\h^ < Ce^"^*. 

The approximated density papp is constructed in order to have the following high order 
approximation: 

dtPapp ~l~ M Papp + Efj^pp.S/ xPapp — Rapp (6-5) 

where Eapp = VA^^{pl^pp + Papp) and Rapp is a remainder satisfying the estimate: 
\\Rapp\\Hi^-2N-i < C(5^+^exp((iV + l)ReAt). 

Let A'^ G N* to be chosen later and take any S > such that S > 2N + 1 . We choose 
also 9 < 1 (to be fixed later) such that ^ > and define ts such that 6 = (5exp (Re Xtg). 

Now we can construct the uj = i j by recursion; we will ensure that for all 1 < j < 

N, J (n+ + uj)dx = and 

< Cexp(jReAt). 

Suppose we have uj for j < k. Then we define Uk+i as the solution of the linear 
equation: 

k 

+ Ei^,_j.VxU,=0 (6.6) 

3=1 

with Ej = VA~^(n^ + uj) and n^.+i(0,x) = as initial condition. Intuitively, n^+i is 
chosen in order to counterbalance the non-linear interaction between the previous terms of 
the expansion. 
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Thanks to Corollary 14.11 with 7 G]ReA,2ReA[, we get the following estimate: 
i-t k 

\uk+i\\HS-{k+i) < / \\e^'^'^*'''\^Ef .VxUk+i-j + E^_^-^_j.\7xUj)\\^L-(k+i)ds 

< C e'^^*~^^(y^ \\Ef\\fj3-{k+i)\\uk+i-j\\HS-k + \\E^_^i_j\\fjs-{k+i)\\uj\\fjs~k)ds 

Jo 

< C [ e^(*-") exp ((A; + 1) Re Xs)ds 

Jo 

< Cexp((/c + l)ReAt) 

Note also that since ^ / (u^^^^^i + u^^^)dx = 0, we clearly have 



Now we can see that: 

dtPapp ~l" ^ Papp ~l~ ^app'^ xPapp — Rapp (^•'^) 

with Rapp = ^2N>j+j'>N^''^'' Ej-.VxUj'. Then, noticing that for all t < ts: 

(5 exp (Re At) < 9 < I, 

the following estimate follows: 

\\Rapp\\HS-2N-i < Cjv(5^+i exp {{N + 1) Re At) (6.8) 

where Cat is a constant depending only on A^. 

Now we consider the solution p to (13. 4p such that p(0) — p, = PappiO)', the equation 
satisfied by w = p — p — papp is the following: 

dtw - + E^.V^w + E^pp.V^w + E^.V^papp = -E^.V^p - Rapp (6.9) 

with Eyj = VA-i(ty+ -W-). 

Then we want to estimate ||w||j;,2 by using some modulated energy inequality. To this 
end, we multiply by w and integrate with respect to x: 

^1^11^2 < J \E:^.VxPappw\dx + J \E:^.V^pw\dx + \\Rapp\\L2\\w\\L2 

< (llV^PapplU- + l|Vx/i||L°°)||-E^il|L2|kl|L2 + ^\\w\\h + ^W^appWh 

< C ((1 + ||V^Papp||L°°)||w^|li2 + \\Rapp\\'i2) 

< C ((1 + \\V.,papp\\L^)\\w\\l2 + Cn6^^^-^^^ exp (2(iV + 1) Re At)) 

But for a > such that 2 + a < 5 — and for t < tg, we can control the Lipschitz 
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norm of papp. 



N N 

< y^(5'iiniiig2+c» < i^'ii^iiijj-s- 



1=1 1=1 

N 



< ^(5^exp(iReAi) 

i=l 

- ^ -1-0-2 



i=l 



Now choose such that: 



3C 

iV+l>— (6.10) 



By Gronwall's lemma we consequently get: 

Hp - M - PappWi^ = \\w\\l2 < CnS^+^ exp {{N + 1) Re Xt) < Cn0^^^ (6.11) 
On the other hand we have a bound from below for the norm of Papp, for t = ts'- 



N 



L2 > (511^111^2 - ^5*11^^11^2 



i=2 

N 

> 5 exp (Re Xts) — (5' exp (i Re Xts) 

1=2 

N 
j=2 

> 

- 2 



Finally we have, for t = tg: 

||/5-P||l2 > llPapplliZ - Hp - ^ - /)app||L2 

> h- CnS^~^^ exp ((TV + 1) Re At) > ^0 - CnO^-^^ 

if 6 is chosen small enough with respect to N. This proves the expected instability result 

Now, in order to prove the exponential growth of the electric field, we use the conser- 
vation of the energy proved in Theorem 15.11 We have for any t > 0: 

\\p{t) - p\\l - J iVVit^dx < \\p{0) - plli. - I iVVm'dx 

(6.12) 

which implies that: 

r??<||p(t5)-//||i2 < J \VVits)\'dx + ||p(0) - plli. 
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We can consider that 6 < 5q < r/i/2, so that: 

I \VVits)\''dx > L{T+ - T-){ril - Sl) := r?i (6.13) 

This proves (j6.2p . □ 

By Remark 15.41 we can be a httle more precise on the growth of the densities. 

Remark 6.3. With the same notations as in the previous theorem, there exists r/s > 
such that: 

\\p^it5)-pi''\\L'^>m (6.14) 

and 

\\p~{ts) - /i"||L2 > r/3 (6.15) 
which means that both the hot and the cold plasma are unstable. 

Proof. According to Remark 15. 4| we have 

\\p~{ts) - /x~||i2 - \\p^{ts) - p^\\l2 = \\p'{0) - P~\\l2 - \\p^{0) - p^Wl^ (6.16) 

which implies that 

\\\p'{ts) - p'\\l2 - Wp^its) - P^Wl^I < So 
By Theorem 16. H we have 

\\p~iU) - P~\\l2 + \\p~^{t5) - P^\\l2 > Vi 
Assuming as in the previous proof that 5q < ei/2, this clearly implies that 

\\p'{ts) - p'\\l2, Wp'^its) - /-i^llia > m/4 := V3, 
which proves our claim. □ 

7 Conclusion 

We have finally managed to provide a mathematical explanation of stability in the "good 
curvature" region and instability in the "bad curvature" region with our simplified nonlinear 
model. In our analysis we have pointed out that large temperature gradients brought 
nonlinear stability even in the bad curvature region. In other terms, if there is enough 
heating, there is good confinement: this is the H-mode. 

A first natural extension to this work would be to generalize the stability/instability 
result to the kinetic model (|3.ip : 

dtf-^\v\^d,J + E^.VJ = 
E = -V^V 

-A,y = J fdv-1 

This shall be the object of a future work. Another important issue is to understand the 
influence of the gyroaverage operator, that we have neglected in this work. 

In "real" tokamaks, the next step towards confinement consists in considering a mag- 
netic field with a variable direction, i.e. B = Bqc^ + Bieg. At leading order, particles 
still follow the magnetic field lines: consequently, with such a twisting field, particles from 
the "bad curvature" region travel every now and then to the "good curvature" region. We 
accordingly expect overall confinement for the plasma. 
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A very challenging and interesting problem would be to prove overall confinement with 
such a twisting magnetic field. In this case, the accurate parameter to consider is the 
so-called safety factor, which stands for the number of times the magnetic field lines twist 
around the torus the long way for each time they twist around the short way. In "real" 
tokamaks, it has to be chosen with precaution in order to get good confinement properties 
(see [27]). But it seems to be much more complicated, since one has to deal with many 
drifts due to the geometry of the magnetic field. 

Finally let us conclude by mentioning that the analysis of confinement provided in this 
paper is very naive since it is now well known that there is a loss of confinement in toka- 
mak plasmas, referred to as anomalous transport. Many models have been proposed and 
intensively studied to justify these phenomena: it would be very interesting to study some 
of them from the mathematical point of view. 
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